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Abstract. BATSE, Ulysses, and TGRS and KONUS on WIND detected four gamma- 
■ ray events within 1.8 days during 1996 October 27-29, consistent with coming from the 

same location on the sky. We assess the evidence that these events may be due to a 
series of bursts from a single source by calculating the probability that such a clustering 
in position and in time of occurrence might happen by chance. The calculation of this 
probability is afflicted by the usual problem of a posteriori statistics. We introduce a 
clustering statistic, which is formed from the "minimum circle radius" (i.e. the radius 
(f) ■ of the smallest circle that just encloses the positions of all the events) and the minimum 

time lapse (i.e. the time elapsed between the first and last event). We also introduce a 
second clustering statistic, which is formed from the "cluster likelihood function" and 
the minimum time lapse. We show that the use of these statistics largely eliminates 
• the "a posteriori" nature of the problem. The two statistics yield significances of the 

clustering of 3.3 x 10 4 and 3.1 x 10 5 , respectively, if we interpret the four events 
as four bursts, whereas the clustering is not significant if we interpret the four events 
as only three bursts. However, in the latter case one of the bursts is the longest ever 
observed by BATSE. 
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>< '• INTRODUCTION 



The question of whether gamma-ray burst sources repeat on timescales of a year 
or less has occasioned much debate since the public release of the first BATSE 
(IB) catalog [1]. Information concerning such repetition would prove enormously 
constraining on gamma-ray burst models. Most models of burst sources at cosmo- 
logical distances cannot plausibly repeat on such timescales, since in these models 
the burst destroys the source. Models of burst sources located in a galactic halo, 
on the other hand, are expected to repeat on such timescales on energetic grounds 
[2]. Thus, confirmation of such repetitions would create many difficulties for cos- 
mological models, while ruling them out would seriously constrain galactic halo 
ones. 

There have been several claims and counterclaims made in the literature on 
the subject. Global evidence for repetition was found in the IB catalog but not 
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FIGURE 1. Sky locations of the four events. The crosses represent the BATSE locations. The 
errors circle radii are computed using the best-fit model for the BATSE location systematic error 
from Graziani & Lamb ( [5]). The IPN location is also shown. 

confirmed in the 2B and 3B catalogs; individual clusters suggestive of repetition 
were also found in the IB catalog, although their statistical significance was not 
impressive. For details, see the discussion in [3]. 

This, then, was the situation on 1996 October 27, when the first of a series of 
four BATSE GRB triggers occurred. By the time of the fourth trigger, 1.8 days 
later, the four burst triggers were conspicuous for having occurred with a relatively 
short time and for having locations consistent with a single joint source on the sky 
[4]. The locations are shown in Fig. 1. The positional and temporal coincidence 
of these events is striking, and provokes the question of whether it is probable that 
they were produced by a single, repeating gamma-ray burst source. We address 
precisely this question in what follows. 

THE DATA 

The four burst triggers (BATSE trigger numbers 5646, 5647, 5648, and 5649) 
occured in two pairs. The first pair occurred about 18 minutes apart on October 
27, while the second pair occurred 11.2 minutes apart 1.8 days later, on October 29. 
The last event was considerably more intense and fluent than the first three, and 
was also detected by KONUS-WIND and by ULYSSES, resulting in the truncated 
IPN annulus shown in Fig. 1. 

While it is not unnatural to regard these four triggers as four gamma-ray bursts, 
the BATSE team has raised the possibility that the last two triggers may in fact 
correspond to one, extremely long burst, in which case there would only be three 
bursts [4]. This possibility is difficult to assess, as there is no model-independent 
way of making the distinction between a gamma-ray burst trigger and a gamma-ray 
burst. At a minimum, the claim that triggers 5648 and 5649 are due to the same 



gamma-ray burst is open to debate. Certainly, the triggers were initially classified 
as independent bursts. There does not appear to have been any emission in the 
approximately 11 minutes between the end of the third event and the onset of the 
fourth, and if the two events were part of the same burst, it would be the longest 
burst observed to date by BATSE. 

Nonetheless, the possibility exists that there were in fact only three bursts ob- 
served during the period in question, and we will consider this possibility in our 
analysis. In the following, when we do consider the case of only three events, we will 
do so by dropping the third event (trigger number 5648), since trigger number 5649 
was much more fluent, so that its location is much more accurately determined. 

ANALYSIS 

We adopt the frequentist statistical methodology of hypothesis testing. The 
null hypothesis under scrutiny is that the four (or three) bursts in question all had 
distinct sources. It is necessary then to identify a test statistic that allows us to ask 
some form of the question "assuming the null hypothesis, what is the probability 
that the observed position-time clustering should have occurred?" 

The chosen test statistic should certainly be sensitive to the kind of clustering of 
interest. However, there exists a danger that the statistic, being chosen a posteriori, 
will be unfairly selected to confer an excessive significance upon this particular data 
set. As an example of this danger, consider the following calculation. The four 
BATSE point locations may be inscribed in a circle of radius 6.1°, and their times 
of occurrence in a time window of duration 1.8 days. Using a rate of BATSE burst 
triggers of 0.8 day -1 , this works out to an expected number of events in our "time- 
solid-angle" window of 4.1 x 10~ 3 . The resulting Poisson probability for seeing 
four events is 1.1 x 1CT 11 . Multiplying by the number of such "time-solid-angle" 
windows in the BATSE catalog up to 1996 October 29, we obtain the apparently 
impressive significance of 6.9 x 10~ 6 . 

The problem with the above calculation is that it would produce greatly different 
significances for different choices of the "time-solid-angle" window. On the one 
hand, there is no unique way to pick the window for this kind of calculation. On 
the other hand, the choice made above is the one that (unfairly) maximizes the 
ostensible significance of the result for the particular data that was observed. 

Guided by the above example, we have sought out time-solid-angle clustering 
statistics that are not arbitrary, in the sense that they do not contain parameters 
which may be adjusted independently of the data to "optimize" the significance 
(such as the radius and duration of the window in the example). We have identified 
two suitable test statistics: the minimum circle radius and the cluster likelihood. 

The minimum circle radius approach is a twist on the previous example. Instead 
of calculating the probability that a 1.8 day-6.1° "time-solid-angle" window in the 
BATSE catalog should contain four or more bursts, we calculate the probability 
that there should be four bursts in the BATSE catalog characterized by a minimum 



circle of radius 6.1° or less, and a first-to-last duration of 1.8 days or less. By 
minimum circle, we mean the unique smallest circle that just includes all of the 
events. By this modification, we eliminate the arbitrariness in the choice of the 
window parameters, which are set by the data itself. The probability may no 
longer be calculated by the naive Poisson method described in the example. It is 
necessary to calculate the distribution function for minimum circle radii for Poisson 
point processes on the sphere. 

We have calculated this distribution function. Suppose the Poisson process has 
a rate R, so that in a time r it averages Rr events. The probability that it should 
produce n events inscribed by a minimum circle of radius 9 or less is given by 
Q n (0o', Rr), where 

Q n (9 ; Rr) = e~« +[Rrn-(n- 1)(„ + 2)] £ ||, (1) 

where n = Rt(1 — cos6* )/2 is the expected number of events in the circle. 

We actually need the F n (6o,To), the probability that n events should inscribe 
themselves in a minimum circle of radius 9 or less and should have a first-to-last 
duration r or less. We have rigorously derived the expression for F. Assuming 
that r <C T, where T is the lifetime of the experiment, we find 

T 

F n (0 , r ) « Q n (8 ; Rt ) x — x n. (2) 

Applying Eq. (2) to the data, we find that assuming there are four events in 
the set, the probability is 3.2 x 10~ 4 , while assuming three events, the probability 
is 1.2 x 1CT 1 . We have confirmed these results by Monte Carlo simulations - 
3.3 x 10~ 4 and 1.1 x KT 1 , respectively. 

We have as yet made no use of the information borne by the BATSE error 
circles. They may be introduced by means of the cluster likelihood, which was 
introduced in a Bayesian framework by Graziani & Lamb [6]. Briefly, it is the 
probability that n events should have their measured locations assuming they have 
a common source, integrated over their unknown source position. In the Gaussian 
error approximation, the expression for the cluster likelihood is 
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where is the location of the ith event, is its Gaussian variance, T 2 = 

E?=i^- 2 andz = r 2 Er =1 ^M 2 - 

The usefulness of C for clustering studies is apparent from its form. It is a product 
of two competing factors: The first, rational factor increases with increasing n 
(recall that the a are generally much less than 1), and so rewards larger putative 
clusters. At the same time, the second, exponential factor penalizes clusters that 



TABLE 1. Summary of Statistical Significances 





minimum Circle Probability 


Cluster Likelihood Probability 


Three Events 


1.2X1CT 1 


1.6xlCT 2 


Four Events 


3.3xlCT 4 


3.1xlCT 5 



are too dispersed. This tension makes C an excellent clustering statistic. Sadly, its 
distribution function is not known, and must be obtained by Monte Carlo. 

We have performed the simulations, creating BATSE catalogs with uniformly 
distributed, non-repeating bursts whose location errors are randomly drawn from 
the BATSE location error distribution. In the case of n events in a cluster, we 
record the number of times that n or more events occur within 1.8 days or less 
and have a value of £ exceeding that of the data. In the case of four events, the 
probability that we obtain is 3.1 x 1CT 5 , while in the three-event case the probability 
is 1.6 x 10" 2 . 

CONCLUSIONS 

Our results are summarized in Table 1. We find that the evidence in favor of 
repetition of burst sources depends very strongly upon the assumed number of 
events. Assuming four events, the evidence is quite strong. In the case of three 
events, the significance is much weaker, but in this case one must accept that one 
of the events was the longest gamma-ray burst ever observed by BATSE. 

We acknowledge support from NASA grants NAGW-4690, NAG 5-1454, and 
NAG 5-4406. 
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